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Abstract 

m . 

^ ^ We study the phenomenology of a supersynimetric scenario where the next-to- 

hghtest superparticle is the hghter stau and is long-hved due to a very weakly coupled 
OO I lightest superparticle, such as the gravitino. We investigate the LHC sensitivity and 

^*| ■ its dependence on the superparticle spectrum with an emphasis on strong produc- 

tion and decay. We do not assume any high-scale model for SUSY breaking but 
on . work along the lines of simplified models. Devising cuts that yield a large detection 

efficiency in the whole parameter space, we determine the LHC's discovery and ex- 
clusion potential. This allows us to derive robust limits on m^^ , 
and mj . We briefly discuss the prospects for observing stopped staus. 
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1 Introduction 

In supersymmetric scenarios where the gravitino is the lightest superparticle (LSP), the 
decay of the next-to-LSP (NLSP) is suppressed by the very weak gravitino interactions, 
as long as R parity is conserved. For a gravitino mass above a few keV, NLSPs produced 
at a collider usually have a decay length that is large compared to the size of the detector. 
The same is possible in axino LSP scenarios. For a charged NLSP such as a stau, this 
leads to the exciting possibility of finding SUSY in events with no missing energy and 
two charged tracks leaving the detector. 

In a collider experiment such a long-lived stau can either be produced directly or in 
a decay chain following the initial production of a pair of superparticles. At the LHC 
the production of the strongly interacting squarks and gluinos clearly has the biggest 
potential to dominate over other production mechanisms. This is why in this work 
we concentrate on the strong production. Unfortunately, a considerable part of the 
more than 100 free parameters of the MSSM influences the signature of staus via the 
appearance of intermediate sparticles in the cascades. This general problem of the large 
SUSY parameter space has often been tackled by studying constrained models such as 
the CMSSM. The nonobservation of SUSY has severely reduced the constrained model 
parameter space. Among other things, this has driven the interest in model-independent 
studies, including regions in parameter space that are not covered by constrained models. 



There are basically two ways of going beyond constrained models. The first is to 
use a well-motivated set of 19 free parameters in the framework of the phenomenological 
MSSM [1] and perform a Monte Carlo scan over the (still vast) parameter space, display- 
ing the behavior of observables in a scatter plot (see e.g. [2]). The second is to reduce 
the parameter space in a bottom-up approach more drastically after identifying the most 
important low-scale parameters that determine the signature. This latter approach finds 
its realization in the so-called simplified models |3lll]- So far the idea of simplified models 
has not been considered in the case of a very weakly interacting LSP. Nonetheless, it is 
especially suitable in the case of a long-lived stau scenario, as we will show in this paper. 

We will determine the LHC sensitivity for a general long-lived staiu scenario utilizing 
a simplified-model approach. We consider the 8TeV and the 14TeV LHC runs. If not 
stated otherwise, we refer to the 8 TeV run. We consider the particles of the MSSM as 
the only particles involved in the interactions inside the collider. We shall assume that 
there is no accidental phase space suppression that renders any sparticle other than the 
stau NLSP long-lived. 

We will define a set of simplified models in section |2] We will divide the problem into 
two parts — the production and the decay. We will first discuss how the production cross 
section depends on the mass pattern in the squark and gluino sector. This discussion is 
general and applies to any LSP scenario. However, we will consider a range of sparticle 
masses around the LHC limits that are typically reachable in the long-lived stau scenario. 
Our results will be presented for the two limiting cases of a common squark mass and of 
a scenario in which one of the stops is much lighter than all other squarks. In each case 
we will allow for two free parameters: the gluino mass and either the common squark 
mass or the stop mass. We will then consider the decay of colored sparticles into the 
stau NLSP and discuss the impact of the intermediate sparticles in the decay chain on 
the observables. This discussion is based on the observation that the direct signature 
of the stau (rather than SM particles from the cascade) provides the most significant 
contribution to a potential discovery or an exclusion. Since the identifiability of the stau 
depends strongly on its velocity, our considerations are driven by the quest of finding 
the quantities the stau velocity dominantly depends on. This necessitates introducing 
the stau mass as the third free parameter. We will define three simplified models that 
serve as limiting cases and thus are able to capture the phenomenology of any realistic 
spectrum within the long-lived stau scenario. 

In section [3] we will study the relevant background sources and discuss the rejection 
obtained by a cut on the velocity of a stau candidate. Due to two distinct measurements 
of the velocity and the requirement of two staus per event, an excellent background 
rejection is possible. For certain regions in parameter space a large number of staus are 
likely to be missed by the current trigger settings. Hence, we will propose a dedicated 
trigger in order to be able to record the corresponding events. In section 13.41 we will 
introduce selection criteria that achieve high efficiencies in the whole parameter space. 

In section |4] we will present the results obtained in a Monte Carlo study by scanning 
over the defined parameter space. We will show the discovery and exclusion reach in 



^The analysis is virtually identical for any charged slepton NLSP. 



terms of the stau, squark and gluino masses, both with and without utiUzing the proposed 
trigger setup. These results represent conservative, model-independent bounds and thus 
can be applied to all models with a long-lived charged slepton. 

Some staus may be stopped inside the detector and decay much later. We will discuss 
the potential for observing such decays in section [5l estimating the number of staus that 
are available for the analysis. 

2 Simplified models for long-lived stau scenarios 

The task of simplified models is to reduce the huge SUSY parameter space. The same is 
true for constrained models, but those obtain a reduction by imposing boundary condi- 
tions with few parameters at a very high energy scale. In simplified models the reduction 
is driven by the signatures at colliders. We aim to determine the most important low- 
energy parameters governing the LHC sensitivity to long-lived staus originating from 
cascades following the strong production of SUSY particles. 

2.1 Production 

Obviously, the masses of the produced squarks and gluinos play an important role — the 
production cross section depends on them. The first simplifying assumption we discuss is 
a common mass for all squarks. In the high-mass region, which we are primarily interested 
in, the production of squarks requires partons with a large momentum fraction. The 
corresponding parton distribution functions (PDFs) of u and d quarks are much larger 
than those of heavier quarks and of antiquarks. Hence, governed by the large contribution 
of the first-generation squarks, ^^ production dominates over qq {see upper right panel of 
figure [8| , and gq dominates over gq. The latter is even negligible and thus not considered 
here at all. 

The effects of abandoning the assumption of a common squark mass are displayed 
in figure [1] The left panel shows the next-to-leading order (NLO) cross section for qq 
production computed by Prospino 2 |5] as a function of ?TT-T(red dashed curve originating 
at 1 on the ordinate) , rriu (green solid curve) and m~ t = m-r = rriu (blue dotted curve) . 
In each case, all other squark masses are degenerate at nig = 1600 GeV. We normalized 
the curves to the qq production cross section at the point where all squarks are mass 
degenerate. The ratio m-g/mq is chosen to be 1.68. This is the value where ^^ production 
contributes maximally to the total cross section of colored sparticles in figure (8] From 
figure [1] it is obvious that the contribution of the first-generation squarks is dominant — 
in the case where d and u masses increase simultaneously, the qq cross section drops 
drastically. Decoupling either only u or only d has a much less drastic effect. Decoupling 
the squarks other than u and d does not significantly change the production cross section. 
Additionally, we plot the analogous curves for ^Q^ production (curves originating at slightly 
above 0.1 on the ordinate). Here, the i-channel contribution, which introduces the flavor 
dependence, is much less important. Consequently, if we decouple both d and u, the 
qq cross section drops much more moderately than the one for qq and hence becomes 
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Figure 1: NLO cross sections for the production of colored sparticles at the 8TeV LHC. Left: Squark- 
squark (curves originating at 1 on the ordinate) and squark-antisquark production (curves originating 
below) as a function of the mass irn of different squark species. All other squark masses are degenerate 
at TTiq = 1600 GeV and m-g/rriq ~ 1.68. The curves are normalized to the squark-squark production 
cross section for the fully degenerate case, i.e., rrii = mq. Right: Total cross section for all production 
channels (normalized to the cross section with degenerate masses) along a typical LHC sensitivity limit 
(see figure [8]) parametrized by the ratio m-g/mq. We decoupled d, u, and both flavors, obtaining the 
red dashed, green solid, and blue dotted lines, respectively. The purple dot-dot-dashed curve shows the 
cross section for q = t\ and all other squarks decoupled. 



dominant. As we stated above, 'gq is negligible compared to the other contributions, 
and so is gq in the case of decoupled u and d — the cross section drops by two orders of 
magnitude. Consequently, for decoupled first-generation squarks, qq and gg remain the 
most important production channels. The latter channel is evidently much less sensitive 
to the squark masses. 

In the right panel of figure [1] we plot the total cross section for all production channels 
(55) IJQi QQ: qq) along an expected 8TeV LHC sensitivity limit (black dashed line in the 
left panel of figure [8)) parametrized by the ratio mg/rriq. As before, different squark 
species are decoupled while the remaining squarks have a common mass rrig, and for 
each value of m-g/mq the cross section is normalized to the cross section with degenerate 
squark masses. The vertical dotted line marks the ratio m-g/vriq ~ 1.68 used in the left 
panel. Stop production is not included in the curves corresponding to decoupling d or u. 
However, it would change only the blue dotted curve noticeably, giving an enhancement 
by a factor 4/3. In each case considered so far, the right- and left-handed squarks were 
treated uniformly. Since the strong interaction is not chiral, decoupling only ^L or ^R 
merely results in a combinatorial factor that does not involve information from the PDFs. 

The lowermost curve in the right panel of figure [1] corresponds to decoupling all 
squarks but the lighter stop. Decoupling all squarks except either ^2 or one 6, c or s 
squark would give the same result. This scenario provides the limiting case of a minimal 
cross section. 

We are left with two limiting setups, the one with a common squark mass and the 
light stop scenario. Each contains two free parameters, nig and ?n,q, or rrig and m^ , 
respectively. We will focus on the former case but come back to the latter scenario in 
section 14.31 From the discussion in this section the LHC reach for a general scenario can 



be estimated. 

2.2 Decay 

After the production of squarks and gluinos, these particles decay via a cascade to the 
stau NLSP. The mass difference between the hghtest colored sparticle (LCP) and the 
stau NLSP determines the total phase space available in the cascade. Thus, it strongly 
affects the kinematics of the stau (and of the SM particle radiation). We will therefore 
consider as a third free parameter the stau mass. 

In order to study the impact of the intermediate sparticles in the cascade, we consider 
limiting cases. Following the considerations in [6l [7] and translating them into the stau 
NLSP scenario we find that within the MSSM there are no spectra for which the LCP 
preferably decays via chains containing more than three intermediate sparticles between 
the LCP and NLSP. In large regions of the parameter space, shorter decay chains give 
the dominant contribution. 

We focus on the impact of the mass spectrum on the stau velocity, which is the most 
important quantity affecting the identification of staus. We consider the limit where only 
massless SM particles are produced in the cascade. For a two-body decay the velocity of 
the daughter sparticle i + 1 in the rest frame of the mother sparticle i is 

If the mother sparticle has velocity /3j in the lab frame, the velocity of the daughter 
sparticle in the lab frame reads 



a. 



'j+i 



(l-/3?)(l-/3Si^) (2) 



\ (1 + AA 



a-e)^ 



where G\lx ^^ ^^^ decay angle in the rest frame of the mother sparticle. Assuming a fixed 
mass gap niQ — m„ (= tjilcp ~ "^ri) ^nd considering an (n — l)-step decajo (mo > m\ > 
• • • > fnn) with a uniform probability distribution for ^^-V^ in each decay (i.e., ignoring 
spin correlations), we can compute the mean of the NLSP velocity /3.„ as a function of 
all masses and the LCP velocity /3q, 

/3n =^n(/5o,"ll,---,"^„_l). (3) 

It turns out that /3„ has one minimum (maximum) at the point given by the mass patteri^ 

n — i i_ 

rrii ~ ttIq" rrin (4) 



^Following [4], we refer to a cascade with n intermediate sparticles between the LCP and the NLSP 
as an 'n-step decay'. For instance, the decay q^ ^ ^t is a, 1-step decay. 
^We checked this expUcitly up to n = 3 but expect it to hold for any n. 
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Figure 2: Contours of constant average velocity /3„(/3o,?tii, • • • ,?ti„_i) for mo = 1400 and m„ = 600 
in arbitrary units. (Tlie overall mass scale is irrelevant. For convenience we chose values that could 
be realistic masses in GeV.) Left: /Sj as a function of /3o and m,i. The dashed vertical line denotes the 
extremum according to Q. For /3o < 0.77 it is a minimum. Above this value it turns into a maximum, 
which is not very pronounced, though, as /Jj is nearly independent of mi . Right: P^ as a function of mi 
and 7712 for /3o — 0.6. The central dot denotes the minimum according to ^. The maxima according 
to (O are located in the lower-left, upper-left and upper-right corners. The limiting case of an effective 
1-step decay chain (n — 2) lies on the borders of the contour plot — the upper border {mi = ttiq), the 
left border (7712 ~ 7773) and the diagonal (7771 — 7772). 
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The extrema ([5]) represent the mass-degenerate hmit and correspond (in the approxima- 
tion we are currently working in) to the direct decay of the LCP into the NLSP. The 
result (|4]) is not surprising: it renders all velocities /3j A to be equal — on average each 
decay gives the same contribution to the velocity of the NLSP. In contrast, in (j5)) one 
decay dominates over the others. The result also implies that the extremal values of /3„ 
lie between those of /3„ if m < n. The m-step decays represent a slice in the space of the 
masses in the n-step decays with n — m masses degenerate. This slice clearly does not 
contain the point (|4]). 

Whether (|4]) is a minimum or a maximum depends on /3o- In fact, (j4]) is a maximum 
only if (3q is very close to the speed of light (see the left panel of figure [2)) . The high 
efficiency obtained in the long-lived stau search (see section [33]) pushes the boundaries 
of the LHC sensitivity to high squark and gluino masses. Hence, they will typically be 
produced rather close to threshold and thus /3o is expected to be significantly below 1, 
at least if LCP production dominates. In this case (|4|) is a minimum of /3„. The right 
panel of figure [2] shows the contours in P^ for a 2-step decay as a function of mi and m2- 

We can now formulate appropriate simplified models. We see that the direct decay 
of the LCP into the stau mediated by a nearly mass-degenerate neutralino as well as the 
mass pattern (j4| are reasonable benchmark mass patterns. We will display our results 



for three choices. 

Model A The 'direct decay' via a nearly degenerate neutrahno (mp ~ ttT'ti] 

Model C The mass pattern (|4| for the 3-step decay LCP — t- X2 ~^ ^ ~^ Xi ~^ '^■ 



Model B The mass pattern (|4)) for the 1-step decay LCP — t- x*^ — ^ ''", ?T^^;o = y/mhcpfn^ 



In all cases we will force the respective branching ratios to provide the desired cascades 
(see section |47T] for details). We will assume symmetric decay chains, i.e., the same cas- 
cade for both LCPs produced in an event. We will briefly discuss the issue of asymmetric 
chains in section 14.21 

The model with pattern (|H) and the 2-step decay LCP — )• X2 ~^ X? ~^ ^ turned out 
to lie completely between B and C concerning the LHC sensitivity shown in section 14.21 
By this we implicitly checked also that the appearance of the heavy SM particle radiated 
in the decay of the heavier neutralino does not change the qualitative picture. Threshold 
effects are expected to be small due to the large SUSY masses o 

3 Background estimation and selection criteria 

For the initial production of colored super par tides, each event contains at least two jets, 
two staus and two taus or tau neutrinos. As the identification of taus it not very efficient, 
we do not include them in the signature. Furthermore, we will see that the background 
rejection can already be saturated with staus and jets alone. 

In the detector, long-lived staus show up as muon-like particles, i.e., charged particles 
usually leaving the detector. They can have a velocity f3 significantly below the speed of 
light, which allows one to distinguish them from muons by virtue of a cut on /3. However, 
in some regions in parameter space many staus are produced with a velocity close to 1. 
These regions typically feature spectra with large mass gaps rriLCP ~ m^^. Requiring 
hard jets can alleviate the drop in sensitivity when /5 approaches 1, allowing a relaxation 
of the cut on /3. However, dropping the /3 cut completely will lead to a substantial loss 
of sensitivity due to a dramatic increase of the (then unsuppressed) muon backgrounds 
Thus, we will always require identified staus. If the staus stem from rather compressed 
spectra the jets are expected to be soft. On the other hand, the staus tend to be slow in 
these cases, so identified staus alone suffice for a good sensitivity. 

Since the SUSY particles are always produced in pairs, the largest significance can be 
achieved by requiring two stau candidates in each event. Therefore, we use the following 
signature for SUSY events: 



*An exception could be a dominant decay chain with Higgsinos whose decay produces a heavy Higgs. 
Here, threshold effects can be somewhat more important but are not expected to change the picture 
significantly. 

^A sufficient background rejection might be achievable without a velocity discrimination if a very 
specific signature of SM radiation is considered. However, such a search would introduce a strong model 
dependence, which is against the idea of this work. 



• 2 high-pT) isolated muon-like particles passing a velocity cut and 

• (optionally) 2 liigli-pT jets. 

3.1 Background kinematics 

As background we consider all relevant SM processes providing two (isolated) muons. 
We will examine the behavior under several kinematic cuts in this subsection. In the 
next subsection we will discuss how the background is further reduced by a cut on the 
velocity. 

We consider the Drell-Yan (DY) production of muons i^Zj^ — )• [i[x) and taus i^Zj^ — )• 
tt\ di-boson production (VF'''VF~, W Z and ZZ\ it production, single t production {tW 
plus t6) and associated Wh production, with jets from initial or final state radiation. We 
calculated the cross section for DY production with FEWZ |8] at next-to-NLO (NNLO) 
accuracy. The cross sections for di-boson production j^, tt and single t production 
[TU| [TT| as well as associated Wh production [121 IE] were calculated via MCFM [H] 
at NLO precision. We generated the events with MadEvent 5 |T5]. To regularize the 
coUinear singularity and gain generator efficiency we imposed the generator-level cuts 
p^ > 60GeV (on both leptons) and p^ > 60GeV (required for at least one h quark) in 
the normalization and event generation of the processes DY and VF6, respectively. The 
resulting cross sections are summarized in table [1] 

We performed showering and hadronization with Pythia 6 jl6j . Since we will impose 
a selection criterion requiring two hard jets, the distribution of the two leading jets should 
be reliable up to very high px- Therefore, for the processes DY, W'^W", single t and Wb 
we include up to two additional jets in the matrix element simulation of MadEvent, 
whereas for the processes WZ, ZZ^ and it we consider up to one additional jet — in 
the latter processes (at least) one of the two leading jets is expected to originate from 
the decay of a heavy SM particle. In the case of Wh^ one of the additional jets in the 
matrix element is allowed to be a 6 jet in order to include the VF66 contribution. In order 
to properly match the different contributions to the inclusive sample, which contains 
jets both from showering and from the matrix element, we applied the MLM matching 
procedure |17) and chose xqcut = •p'^ '™™ = 30GeV and QCUT = 40GeV. We used the 
Cteq6l1 PDF set [18]. 

We passed the output of Pythia to the detector simulation Delphes 1.9 [19] and 
applied a cut and count analysis on the Ihco output of Delphes. The reconstruction 
efficiency for each muon was set to 0.9. The trigger efficiency was conservatively set to 
100% for the background. Figure |3] shows the cross sections for all considered background 
processes as functions of various cuts on characteristic variables. These variables are the 
transverse momentum of the muon -p^ and of the jet •p'^ , the difference of the px of the 



®If both of the two hardest jets did not originate from the W or Z, respectively, but from initial state 
radiation, this process would be just another correction to the DY process with an additional suppression 
by the weak coupling due to the production of the extra vector boson. Hence, we can easily estimate 
that such a process cannot compete with the DY process, independent of the cuts we apply on the two 
jets and muons. 
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Figure 3: Inclusive cross section of the considered SM background processes as a function of various 
cuts at the 8 TeV LHC. In all plots we imposed the isolation cut /^ei < 0.2 and the pseudorapidity cut 
\i]\ < 2.4 on each muon. Top left: Requiring at least two muons with pi^ > p^'™'" each, p^'™'" is varied. 
Top right: Requiring at least two muons with p!^ > 80 GeV and Ap!^ > Z\p!^'™'". Z\pJ^'™'" is varied. 
Bottom left: Requiring at least two muons with p^ > 80 GeV and additionally at least two jets with 
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> pi^''""" each. We vary p^''™'". Bottom right: Requiring at least two muons with pi^ > 80 GeV 



and additionally at least two jets with pJ^ > 150 GeV each as well as Ap!^ > Apr. 
Z\pJ^*'""". The key in the bottom right panel holds for all panels. 
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NLO 
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NLO 
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NLO 
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NLO 
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< 0.0002 
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0.0069 
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0.0054 
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4.8 X 10^^ 


8.2 X 10-8 
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Table 1: Cross sections for the considered background sources at the 8TeV LHC. We display the total 
cross sections as computed at the given precision (see text for details) as well as the cross sections times 
efficiency obtained by applying the cuts selection criteria 1-3 defined in section \3A\ but without cuts on 
the velocity /3. The cuts on the velocity belonging to selection criteria 1-3 are only applied in the very 
Icist line. 

two considered muons, 

^P^ = |P^''-P^''|, (6) 

and of the two considered (hardest) jets, 

M'* = K*''-P^?''|. (7) 

To reject the QCD background we require isolated muons, i.e., /rei < I^i'^, where 

E track 1^ v~^ iriCAL 
^rel = Ji \P) 

Pt 

is the relative isolation. In ([S| the sums are performed over all objects within a cone of 
a given AR = y'Arf + AcI)^ = 0.3 around the muon|j 

Except for the Wb + X contribution the background is not sensitive to the precise 
choice of /™f^ in a range from 0.02 to 1. The same is true for the typical signal we will 
consider. In contrast, Wb + X drops significantly with smaller /™f^ and is practically 
irrelevant for the chosen value of /^f^ = 0.2. Hard muons originating from bb are 
suppressed even more strongly by the isolation cut and do not play any role here. The 
same is expected for other QCD backgrounds. 

The DY production is the dominant background providing two hard isolated muons 
even with the requirement of two hard jets. Without the requirement of hard jets, tt+X 
contributes a few percent of the background, while with this requirement its contribution 
constitutes up to approximately 15%. The other sources are small and do not exceed a 
few percent in total. 



^The energy deposition and track associated with this particle itself are excluded from the sums. 
Throughout this work we use the anti-fct jet clustering algorithm |20| . 
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3.2 Discrimination via the velocity 

The main difference between staus and muons is tlieir velocity, wliicli can be measured 
independently via the ionization loss in the tracker {dE/dx) and via a time-of-ffight 
(ToF) measurement. The relative uncertainty of the ToF measurement of muons (with 
/3 ~ 1) is approximately 0.048 in the ATLAS detector [21] and around 0.06 at CMS 
|22) . The relative uncertainty of the velocity measurement via ionization loss is smaller, 
around 0.035, but it is biased due to truncation effects |23- However, in general it will 
be possible to correct for this bias in an event-by-event analysis. The combination of the 
ToF and ionization loss measurement |22| l2l] yields an unambiguous and very robust 
measurement of the velocity. Hence, we refer to this combination whenever possible. We 
estimated the relative uncertainty of the combined measurement by taking the weighted 
mean of the respective uncertainties for CMS, yielding a^o ~ 0.032. 

A smearing of the velocity of muons is not included in the detector simulation 
Delphes. Assuming that there is no correlation between the velocity mismeasurement 
and the other observableqfj considered in the previous subsection, we treat the back- 
ground rejection due to the velocity cut separately from the application of the kinematic 
cuts discussed there. In other words, for the background we first apply the cuts on /rei, 
'Hi Pt 1 Pt ^Pt ^^"^ ^p^r^ ^° ^^^ generated events and then multiply the resulting cross 
section by the background rejection factor r^ due to the velocity cut. To estimate r^ we 
assume a Gaussian smearing of the velocity with the respective width a'o'^ given above. 
In the case of the signal we refrain from smearing the velocity and use the generator-level 
values (which were passed through Delphes) to allow for an event-based application of 
the cuts. 

Since we always consider two stau candidates, we can combine the velocities of the two 
stau candidates in different ways to formulate appropriate cuts, which lead to different 
factors r^. In the following, we denote the background rejection factor due to a cut 
on a single muon by f^. If the velocities of the two staus within one event tend to 
be correlated, a cut on both stau candidates with the same /3niax will yield the highest 
sensitivity. We will denote this cut as 

/3°</3Sax- (9) 

Since the mismeasurement of the velocity of two background muons in one event is not 
correlated, this yields a background rejection factor of r^ = f'i. If, in contrast, the 
velocity of the two staus in one event is strongly uncorrelated, the cut 



1 - V(l-/3i)2 + (1-/32)2 EE /3° < /3;j,,, (10) 

yields a higher sensitivity on the signal. The background rejection factor is r^ = r^ in 
this case (for the considered Gaussian smearing). 



This is an obvious assumption for the background muons because a deviation from /3 = 1 does not 
originate from the physical process. A dependence of the velocity measurement on the pseudorapidity 
was reported to be small |22) . 
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3.3 Lower limits on the velocity 

For staus with /? < 0.6 the efficiencies of the current triggers at ATLAS and CMS drop 
significantly. In order to improve the trigger for very slow staus the tracker data has to be 
buffered in order to allow for a recording of the tracker data in delay. Very slow staus are 
typically produced in scenarios with large m^^ . We therefore propose a recording of up 
to about four bunch crossings after a trigger by muon-like particles with p^ > 300 GeV 
(500 GeV) at the 8TeV (14TeV) LHC run. We do not expect this to cause the recorded 
event rate to grow significantly. In section 13.41 we will introduce several selection criteria 
one of which assumes the proposed trigger while the others simply require /? > 0.6. 

Although staus suffer a high energy loss due to ionization of the detector material, 
they will often lose only a small fraction of their total energy. Consequently, their ve- 
locity will stay approximately constant when passing the detector. However, since the 
ionization loss increases with decreasing velocity the ionization loss will become relevant 
for the kinematics of the stau if the velocity falls below a critical value — staus might then 
lose their kinetic energy completely and become trapped inside the detector. Since the 
traveling range of charged particles in matter increases linearly with their mass, heavier 
staus are more likely to pass the detector than lighter ones with the same velocity. On the 
other hand, very slow staus typically appear only if they are very heavy, so in conclusion 
stopped staus are rather the exception than the generic scenario^ 

However, to be able to record the tracks of the staus we have to make sure that at 
least one stau reaches the muon chambers to fire the muon trigger!^ Besides, both have 
to pass the tracker to allow for an ionization loss measurement. In order to account 
for this, we used the approximate traveling range of a charged particle in the detector 
material given in |25) (see section |5] for details) and determined the minimal velocity 
that still ensures the required traveling range R as a function of the mass of the stau, 
'^min("^n)- The range R was conservatively set to 1200 gcm~^ and 12 000gcm~^ for the 
requirement of passing the tracker and muon trigger, respectively. For a homogeneous 
detector material with density p = 8gcm~^, this corresponds to a path length of 1.5 m 
and 15 m, respectively. 

3.4 Selection criteria 

We will now introduce a set of three selection criteria, each of which provides strong 
background rejection and high signal efficiency in its domain in the considered parameter 
space. The selection criteria are chosen in a complementary way such that the union of 
these criteria will lead to high efficiencies throughout the whole parameter space of the 
simplified models introduced in section [2] 

In the following, each stau candidate is understood to pass the isolation criterion 
/rel < 0.2 and to lie within a pseudorapidity range |?/| < 2.4. When we cut on the 



®We will examine the possibility of observing stopped staus at the LHC in section [S] 
"This n 
the trigger 



^"This requirement can of course be relaxed if an event contains enough SM particle radiation to fire 
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transverse momentum p^*, we require two particles i that both pass the cut. The values 
given below are valid for the 8TeV (14TeV) LHC analysis. The selection criteria are: 

1. Two stau candidates passing the muon chambers, 

Pt^ > 80GeV (240 GeV) , Api^ > 50GeV (70GeV) , /3° < 0.86 , (11) 

and two jets, 

p'^^"^ > 200 GeV (400 GeV) . (12) 

One stau candidate is required to fire the muon trigger 'in time', 

/5>0.6. (13) 

2. Two stau candidates passing the muon chambers, 

pI^ > 150 GeV (360 GeV) , /3° < 0.88 , (14) 

and 

/? > 0.6 (15) 

for one stau. 

3. Two stau candidates passing the tracker, 

Pt^ > 300 GeV (500 GeV) , p" < 0.73 , (16) 

one of which has to pass the muon trigger chambers, 

/3i > f^'^f^^-^'im^,) , /32 > A^ir (m,J . (17) 

This selection criterion assumes the modified trigger setup proposed in section 13.31 

Figure |4] shows the efficiencies of the selection criteria 1 to 3 for exemplary mass slices 
of the simplified models A, B and C. For all simplified models, selection criterion 1 is 
most efficient if the stau is sufficiently light. The large mass difference between LCP and 
stau ensures the production of high-j5T jets that enable a very good background rejection 
already in combination with the relatively loose velocity cut on /3°, thus cutting away a 
smaller part of the signal than with the cut on /?'-'. For heavier staus, high-pT jets are 
no longer guaranteed. Consequently, selection criterion 2 becomes more efficient, relying 
on the cut on /3^ to discriminate against muons with a mismeasured velocity. If the stau 
is very heavy, many events will contain slow staus that do not pass the cut /? > 0.6. 
However, due to their large mass they have a very large pT- Hence, selection criterion 3 
is optimal for this part of the parameter space. 

The background rejection obtained with these cuts is summarized in table [1] for the 
8TeV LHC. The relative importance of the background sources is similar for the 14TeV 
run. The 14 TeV cuts provide a stronger background suppression as required by the larger 
integrated luminosity considered. 
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Figure 4: Efficiencies of tlie selection criteria 1 to 3 as functions of rrif-^ for fixed squark and gluino masses 
TTig = nig = 1600 GeV (3600 GeV) for the 8TeV (14TeV) LHC and for the simpfified models A {top 
panels), B [middle) and C [bottom). Trigger and reconstruction efficiencies are included. The solid curve 
denotes the union of all three cuts, i.e., the efficiency resulting from selecting all events satisfying at least 
one of the selection criteria. 

4 Exploring the parameter space 

4.1 Computation of the signal 

To estimate the projected LHC reach for the simphfied models described in section |2] we 
performed a full-fledged Monte Carlo study. Here we briefly sketch the computational 
steps of the analysis. The analysis was performed for the 8TeV and 14TeV LHC run. 
For the sake of saving computing time we factorized production and decay. 

For the production we built up a grid of generator-level event files in the m-g-rriq 
plane. We considered the production channels Ijg, gq, qq and qq with a common squark 
mass for q = u,d, s, c, b. We computed the production cross sections for the different 
channels at NLO precision via Prospino 2 and simulated the events with the tree-level 
generator MadEvent 5. Since the size of NLO corrections differs significantly between 
the considered production channels, we performed a channel- wise normalization, i.e., we 
treated each production channel separately throughout the analysis chain. We generated 
a total of 30 000 events per mass point, apportioned between the four production channels 
according to their fraction of the total cross section. We used the CteqGlI PDF set. 

In a second step we passed the MadEvent events to Pythia 6 to perform the decay 
of the SUSY particles (and the showering and hadronization of SM particles). For each 
point on the m^-rriq grid, this allows for a variation of the spectrum below the LCP. We 
computed the decay widths and branching ratios (BR) via SDECAY |26] . The minimal 
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decay chain LCP -^ }^ ^ t was obtained by decoupling all other SUSY particles. The 
decay chain LCP — ?• X2 ~^ ^ ~^ Xi ~^ ''^i ^^s enforced by computing the BR down to the 
X2 via SDECAY and adjusting the following BR accordingly. 

The use of Pythia 6 for the decay of SUSY particles implies certain approximations 
whose validity we have to justify. Pythia 6 factorizes the cascade into decays of on-shell 
particles, using the narrow width approximation (i.e., each decay width is much smaller 
than the corresponding mass difference) , and neglects spin correlations of fermions in the 
chain. Since we are interested in a systematic scan of the free parameters of the simplified 
models including regions where masses of sparticles in the decay chain are nearly mass 
degenerate, we have to ensure the validity of the use of factorization here. 

The only kinematical cuts this analysis strongly relies on concern the (3 of the staus, 
the pt of the staus and the px of the two hardest jets and combinations thereof. Consider 
the case where two SUSY particles in the cascade are very close in mass, i.e., their masses 
are much larger than the mass difference and much larger than the mass of the radiated 
SM particle. In this case the daughter sparticle basically inherits the kinematics of the 
mother sparticle and the SM particle appears as soft radiation. Its distribution might 
not be described well by the simulation. However, we do not cut on SM particles other 
than the jets. As the jets have to survive very severe px cuts, jets from such degenerate 
decays are very unlikely to contribute to our signal. We have explicitly checked the px 
distributions of the two hardest jets and the px and /? distributions of the staus against a 
full matrix element simulation by Whizard |27) and found reasonably good agreement 
for different setups in which we considered the decay q — ?• q^ — t- grri for (nearly) 
degenerate rriq and rriyO , as well as nearly degenerate m^o and m^f^ . The results differ by 
at most 10%. The resulting error on the efficiencies of the selection criteria which enters 
our final results is estimated to be even smaller, and is especially small compared to the 
dominant uncertainties, which are PDF and scale uncertainties that are introduced via 
the production cross section (for an error estimation see section |4^ . Thus, factorization 
is well-justified when relying on the considered observablesO 

We included up to one additional jet in the matrix element. This introduces a source 
of potential double counting. Gluino-gluino production with an additional jet contains, 
for instance, a diagram with an intermediate squark. If the squark is on-shell, this 
contribution is equivalent to the lowest order gluino-squark production followed by the 
decay of the squark. To account for this, Prospino removes contributions from on- 
shell intermediate squarks and gluinos. Accordingly, we removed the same diagrams in 
the event generation in the course of the matching procedure in PYTHIA. As for the 
background, we applied the MLM matching procedure in order to match properly the 
jets from the matrix element and showering. 

We passed the output of PYTHIA to the detector simulation Delphes 1.9. To account 
for long-lived staus we applied minor modifications on Delphes. The reconstruction 
efficiency for each stau was set to 0.9. We assumed a trigger efficiency of 90% |28) . 



^^Relative angular distributions are affected more strongly by the applied approximation. Thus, larger 
deviations are present in AR. However, the dependence of our results on I^d is very weak, as we stated 
in section [H7T1 Consequently, the results are not affected noticeably. 
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4.2 LHC reach for a common squark mass 

We estimate the LHC's sensitivity to observe or exclude the introduced simphfied models 
as a function of the free parameters rrig, rriq and m^f-^^. We consider a common squark 
mass mg in this section. As discussed in section [2?T] the production is dominated by the 
first-generation squarks. Consequently, the derived limits can be interpreted as limits 
on the masses of these squarks. In section W^ we will consider the case of a light stop 
dominating the production cross section. 

As shown in figure |4] the efficiencies of the selection criteria are typically about 0.5. 
On the other hand, as shown in table [U the background expectation is reduced by these 
cuts to less than 10~^ events. Accordingly, a 95% C.L. exclusion^ can always be claimed 
if no events are observed while three are expected. That is, the background rejection 
is saturated with these cuts. In such a situation the analysis is not very sensitive to 
the precise background cross section anymore (see e.g. |29|). Furthermore, discovery can 
also be claimed on the basis of very few events and the results for the exclusion curves 
represent an (typically even conservative) estimation of the 5<t discovery potential. These 
are typical features of the search for heavy stable charged sparticles. 

Figures |5] and |6] show the resulting sensitivity for the 8 TeV and 14 TeV LHC run, 
respectively, for the three simplified models introduced in section 12.21 We visualize the 
variation of rrig^ rrig and rrir^ by showing slices of the parameter space. In the plots 
showing rrig-m^^ and rriq-m^f-^ planes, a fixed ratio rriq/mg is assumed. In the plots of the 
m-g-mq plane, we draw the sensitivity curves by conservatively choosing the stau mass 
that yields the smallest sensitivity at each point of the plane. In addition to strong 
production and decay, we include the production of staus by the direct DY process. In 
order to derive conservative limits we considered the stau mixing angle that yields the 
smallest cross section [S^lifl This contribution is always present and depends on the stau 
mass onlyo 

In model A (blue dashed lines) the decay LCP — t- ^ leads to hard jets and poten- 
tially highly boosted staus. For moderate mass gaps rriLCP ~ "^n the staus are well- 
distinguishable from muons. The additional jet signature leads to a slight enhancement 
of the significance. For larger ?nLCP — ""^ri a large number of staus are rejected by the 
velocity cut and the significance drops sharply, despite the fact that the jets become 
harder. This effect would hide the scenario very effectively from our selection criteria if 
it were not for the direct production, which increases the sensitivity for lower m:^-^ (see 
figure m where we did not include the direct production). Thus, DY production allows 
us to cover the parameter space with a mostly model-independent search. Without it. 



^■^We apply the CLs method here. 

^■^A significantly larger contribution of directly produced staus is possible, see e.g. |30) . 

^*Similarly one could think of including the direct production of the intermediate sparticles that appear 
in the cascades. However, since there is no Drell-Yan production for pairs of pure binos the cross section 
can always be rendered negligible by an appropriate choice of the neutralino mixing. Hence, at least in 
the minimal decay chain, there is no such conservative minimal contribution. Including a minimal direct 
production cross section of the intermediate sparticles in longer decay chains would have a certain effect 
on the results in some of the considered regions in parameter space. However, this would run counter to 
the idea of this work. 
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Figure 5: Projected LHC sensitivity (95% CLs exclusion and approximate 5(7 discovery reach, see text) 
for the models A (blue dashed), B (green solid) and C (red dot-dashed), as well as A (cyan dotted) and 
C (purple dot-dot-dashed) for a reduced set of selection criteria (see text for details). A common squark 
mass rfiq is assumed. In the lower right panel the curves represent the minima in the sensitivity with 
respect to the variation of m^^ . 
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Figure 6: Projected LHC sensitivity (95% CLs exclusion and approximate 5a discovery reach, see text) 
for the models A (blue dashed), B (green solid) and C (red dot-dashed), as well as A (cyan dotted) and 
C (purple dot-dot-dashed) for a reduced set of selection criteria (see text for details). A common squark 
mass rriq is assumed. In the right panel the curves represent the minima in the sensitivity with respect 



to the variation of irir^ . 



we were forced to introduce dedicated searches for each occurring topology to be able 
to cover the small m^-^ region. For stau masses just above the region of dominant DY 
production (and for rrig > rrig) the sensitivity for model A reaches a minimum (see up- 
per right and lower left panel of figure |5|) . Higher luminosities push the corresponding 
mass reach up and cause the production to be closer to threshold. As a consequence, 
the minimum disappears in the projected sensitivity for the 14 TeV 300 fb~ run (see left 
panel of figure [6| . The cyan dotted curves show the sensitivity after dropping selection 
criterion 1, i.e., without taking the jet signature into account. The diminished sensitivity 
illustrates the importance of the additional jet signature in this region. 

Model C has a sensitivity minimum in the intermediate range of mLCP ~ "^-ti ■ Many 
staus are rejected by the upper velocity cut in this region. Furthermore, compared to 
model A fewer hard jets are produced that could compensate for this effect. 

As expected from the discussion in section 12. 2| either ^ or C yields the minimal 
sensitivity for all values of m-g, rrig and rrir-^ . Model B behaves more moderately, resulting 
in a sensitivity in between those of the other models at most points. 

For small mass gaps mLCP ~ "^fi ~^ all three models give the same efficiency. The 
mass pattern of the intermediate sparticles does not play a role in this regime. The staus 
basically inherit the velocity and angular distribution of the produced colored sparticles. 
In this parameter region lower limits on the velocity become important. As stated in 
section 13.31 current trigger restrictions cause the loss of events in which both staus have 
a velocity /3 < 0.6. For small mass gaps this loss is significant. This is illustrated by 
the purple dot-dot-dashed curves that show the sensitivity for model C after dropping 
selection criterion 3, which assumes buffering of the tracker data in order to be able 
to record several bunch crossings in delay. Especially for the 14 TeV run and 300 fb~ 
luminosity, the implementation of such a trigger enhances the sensitivity significantly 
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(see left panel of figure [6| . 

In all the plots, the LHC sensitivities for the simplified models A, B and C span a 
relatively narrow band although the mass patterns of the models are radically different. 
Furthermore, the overall dependence of the sensitivity on m:^-^ is moderate. This situation 
is very different from the one in a missing i?T search like in neutralino LSP scenarios. 
There the sensitivity depends much more on the intermediate spectrum and on the mass 
of the LSP and it is more difficult to cover the limiting cases with appropriate simplified 
models. This shows that the simplified model approach is very suitable for the long-lived 
stau scenario. In the m-g-mq plane, where we took the most conservative choice for m^-^ 
at each point, the simplified models A-C lie even more closely together. From the m^-rriq 
plane plots we can derive conservative projected limits on the gluino and squark masses 
in the common mass scenario. With j C = 16 fb~ at 8 TeV we expect gluino and squark 
masses of m-g < 1.4 TeV and rriq < 1.6 TeV to be either excluded or discovered. With 
jC = 300 fb~ at 14 TeV we are sensitive to m-g < 2.6 TeV and niq < 3.3 TeV. These 
limits allow for a completely model-independent interpretation with respect to all SUSY 
parameters that are not specified in this setup. 

As already mentioned in section [2?2| asymmetric decay chains, i.e., one short and one 
long decay chain in one event, will also appear in realistic models. The selection criteria 
we imposed are not dedicated to such chains and asymmetric cascades will partly fail to be 
selected by these criteria. However, we expect that a dedicated extension of the selection 
criteria will provide equally high signal-to-background ratios. One could for instance 
require one very hard jet and an even stronger asymmetry in the stau kinematics. 

Another aspect we noted in section |2T2] is the presence of heavy SM particle radiation 
in the 2-step decay LCP— )■ X2 ~^ Xi ~^ ^- To check whether threshold effects might affect 
the sensitivity significantly, we computed the corresponding curves for this 2-step decay. 
We explicitly allowed for the 3-body decay of the heavier neutralino into the lightest 
neutralino, which occurs once m~o — rrir-o < niz- The sensitivity curves for the 2-step 
decay lie completely between the 3-step decay (model C) and the 1-step decay (model B) 
in all plots of figures |5] and |6l 

As another check of the generality of the introduced simplified models, we considered 
an inverted mass hierarchy in the simplified model C, mj < ^^yO- In this case a 3-body 
slepton decay occurs. In order to perform these computations we used an extension of 
the SDECAY package |31) allowing for the 3-body decays of sleptons. We found that 
for these spectra the picture drawn here does not change. 

4.3 LHC reach for a light stop 

As a complementary limit we will now consider the case of a light stop ti and all other 
squarks decoupled. We set the gluino mass to m-g = 3m^ . In this setup processes other 
than stop-antistop production are negligible. Although the signature of the decaying 
stops might potentially provide a larger significance with a dedicated selection criterion, 
the benefit is expected to be marginal due to the high efficiencies that are already provided 
by the introduced cuts. Therefore, we refrain from introducing a dedicated selection here. 
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Figure 7: Projected LHC sensitivity (95% CLs exclusion and approximate 5a discovery reach, see sec- 
tion [4]2]) for a light stop fi and decoupled squarks and gluinos otherwise. As in figures [5] and [6] the lines 
denote the models A (blue dashed), B (green solid) and C (red dot-dashed), as well as A (cyan dotted) 
and C (purple dot-dot-dashed) for a reduced set of selection criteria (see section f42l for details). 

Figure [7] shows the sensitivity in the m^ -m^^ plane. According to the lower produc- 
tion cross section for stops, the stop masses that are in reach of the LHC are smaller 
than the accessible squark masses that we have found in the common squark mass sce- 
nario considered in the previous section. The stau mass that is accessible via direct stau 
production is the same, of course. Correspondingly, the gap between the LCP mass and 
the lowest stau masses for which the production via cascades is dominant is considerably 
smaller. Therefore, the effect of the intermediate spectrum becomes less important and 
the curves for models A-C lie even more closely together than in the previous section. In 
particular, the minimum of the sensitivity of model A just above the region of dominant 
direct DY production has disappeared. The velocity distribution of the staus varies less 
significantly with the stau mass in this scenario. 

For mj > m-g we find the same sensitivity to nig as in the large-?7iq limit in the 
previous section. In conclusion, with JC = 16fb~ at 8TeV we expect gluino and stop 
masses of m-x < 1.4 TeV and mr < 950 GeV to be either excluded or discovered. With 
fC = 300 fb~ at 14 TeV we are sensitive to rrig < 2.6 TeV and ttij < 2 TeV. 

4.4 Uncertainties 

We shall briefly discuss the theoretical uncertainties of the cross section computations 
here. We expect these uncertainties to be the most important ones. Further uncertainties 
arise from the event generation and from the simplified detector simulation. We only 
consider the error implied by the scale dependence of the production cross section. The 
error introduced by PDF uncertainties is roughly of the same order. Uncertainties in the 
strong coupling Og are expected to be somewhat smaller. A detailed discussion of the 
errors relevant for the production of colored sparticles at the LHC can be found in |32) . 
In the left panel of figure |8] we exemplarily show again the 8 TeV LHC sensitivity 
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Figure 8: Left: LHC sensitivity in the m-g-mq plane for model B (black dashed curve, same as the green 
solid curve shown in the lower right panel of figure[S]) and its error band according to the renormalization 
and factorization scale variation (green shaded region). Upper right: Fractional contributions of the 
considered production channels to the total production cross section of colored sparticles, parametrized by 
the ratio m-g/mg, which is varied along the sensitivity limit shown in the left panel. Lower right: Relative 
errors in the cross sections due to scale variation for the considered production channels, again changing 
mg/rrig along the sensitivity limit. 

for model B and a common squark mass (black dashed line). The green band shows 
the uncertainty implied by varying the scale in the range m/2 < fi < 2m, where ^ = 
^^F = fJ-R is the factorization and renormalization scale and m is the averaged mass of 
the produced sparticles. The error from scale dependence translates into uncertainties 
of roughly ±20 GeV to ±40 GeV in the squark and gluino masses. The lower right panel 
of figure [8] shows the relative error |(Tm/2 ~ o'2m|/2o"m of each production channel as a 
function of the ratio mg/rrig, which is varied along the sensitivity limit shown in the left 
panel. The upper right panel of figure |8] shows the ratio of the respective cross section 
to the total cross section, again along the sensitivity curve. The average relative error 
is around 0.2 and is lowest in the region of a dominant qq channel. The gg production 
involves very large uncertainties for large rrig/mq. However, in this region the cross 
section for gg production is very small. 



5 Stopped staus 

So far we considered the signal of staus as heavy stable charged particles passing the 
detectors. As we already mentioned in section 13.31 very slow staus can lose their kinetic 
energy completely and get trapped inside the detector. These staus will then decay into 
the LSP inside the detector. This causes a signal that can be recorded with dedicated 
trigger algorithms |33 | I34 | [35]. The decay of the stau reveals very interesting information 
about the LSP and possibly even about the origin of SUSY breaking. Once the mass of 
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the stau is known ^i one can determine the LSP mass from reconstructed 2-body decays. 
Furthermore, for a gravitino LSP an unequivocal prediction of supergravity can be tested, 
namely the proportionality of the stau lifetime 

ASirml ,„ Mp / mLr, \ ~ 
m~ \ m~ I 

Ti \ ri / 

to the Planck mass squared |36| l37] . In scenarios with an axino LSP, stau decays may 
provide insights into the Peccei-Quinn sector |38) . 

In this section we will estimate the prospects for the LHC to observe decays of stopped 
staus by computing the number of staus that are expected to be stopped inside the 
detectors in the framework of our simplified models]^ The mean range i? of a charged 
particle (defined as the average thickness of absorber material the particle is capable of 
traversing) grows linearly with the particle's mass. Its dependence on 7/3 is nearly linear 
in the double logarithmic plot |25) in the region of interest 7/3 < 1. In fact, 

^°^io ( ^I/gTv' ) = '^i + "2 logio (7/3) (19) 

approximates R{'yP) to a precision better than 1%. For iron ci ~ 2.16 and C2 — 3.32. 
From (|19p we determine the maximal velocity of a stau with mass m^^ to be expected 
to stop inside the detector, /3max("^ri)- Conservatively, we assume a maximal range of 
-Rmax = 2400 gcm~^. This is a conservative estimation of the thickness of a CMS-like 
detector in the central regionj^J 

Figure [9] shows the expected number of events that provide at least one stau which 
is expected to stop inside the detector, i.e., with /3 < /3max("^Ti)- According to the 
discussion in section 12. 2| for large mass gaps rriLCP ~ "^n staus are expected to be 
considerably slower in model C than in model A. As a consequence, the numbers of 
stopped staus typically differ by an order of magnitude. 

At the 8TeV LHC only a few stopping events are expected in the mass region of 
interest. Thus, most likely it will not be possible to study the properties of stau decays 
in detail. The high-luminosity 14TeV run, however, will provide reasonable numbers of 
stopped staus in scenarios that lie within the discovery reach of the full 8 TeV dataset or 
the early 14 TeV dataset. Consequently, such scenarios are expected to be accessible in 
sufficient detail to determine the stau lifetime. Determining the LSP mass seems feasible 
as well, unless m-LSP ^ 0.1 m^-^, in which case the mass determination becomes extremely 
difficult because it requires a very precise measurement of the tau recoil energy. 



^^The mass is determined already at the stage of discovery by measuring the momentum and the 
velocity [211 [22] . 

^^Since a detailed study of the LHC's potential for measuring the stau lifetime and LSP mass requires 
a precise detector simulation taking into account details of the LHC operating schedule, we leave this to 
the experimental collaborations. Combining our results with those of |34) it may be possible to estimate 
the prospects for measuring the lifetime in the simplified models. 

^^We took into account 20 layers of silicon (2.33 g/cm^), each 0.5 mm thick; 1 layer of ECAL crystal 
(8.28 g/cm^), 23 cm thick; 16 layers of HCAL brass (8.53 g/cm^), each around 6cm thick; 1 magnet of 
NbTi (5.6g/cm^), 31.2 cm thick; 3 layers of iron yokes (7.87 g/cm^), 30, 63 and 63 cm thick [35]. 
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Figure 9: Number of events that contain at least one stau which is stopped inside the detector for 
the simplified models A (blue dashed), B (green solid) and C (red dot-dashed). The iso-event-number 
curves are displayed in the m-g-mr-^ plane where we have chosen the common squark mass scenario and 



A more detailed study of stau decays including the measurement of the spin of the 
LSP via reconstructed 3-body decays |36| l38] may only be possible at an e'^e~ collider 
with a dedicated detector environment, see e.g. |40| W\\ W2\ 143) . 



6 Conclusions 

We have studied the potential of the LHC to discover or exclude scenarios with a long- 
lived stau (or another charged slepton) in a simplified-model approach. The production 
of colored sparticles is likely to be the dominant production mode at the LHC. SUSY 
events are characterized by muon-like particles that leave the detector and can be slow. 
In addition, hard jets can be expected in some cases. We have defined 2x3 discrete 
simplified models covering the limiting cases concerning the production and the decay, 
respectively. Each model contains only three free parameters, m^^^ rrig and either the 
common squark mass nriq or m^ . We have also included the Drell-Yan production of stau 
pairs. 

Due to the very prominent signature of long-lived staus, exclusion and discovery take 
place on the basis of a very few events. Consequently, a discovery could be established 
in a rather short period of time. Due to the direct DY contribution, regions where both 
staus are typically too fast to be identified (where SM particle radiation would provide 
a more significant signal) are not present. For the same reason, no very specific cuts are 
required, enabling a model-independent analysis. In other words, it is possible to cover 
the whole parameter space with a small number of selection criteria that yield both a high 
signal efficiency and a very good background rejection. In almost the whole parameter 
space, we have found a signal efficiency around 50%. Even in the most challenging region, 
the efficiency does not drop below roughly 20%. This shows that in the long-lived stau 
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scenario no regions exist where the theory effectively hides from observation. This is a 
striking difference to the neutrahno LSP scenario, where both compressed spectra and 
very stretched ones are very hard to observe. 

If all squarks have a common mass the most conservative projected limits for JC = 
16 fb-^ at the 8TeV LHC (/£ = SOOfb"^ at 14TeV) are rUg > l.GTeV (m, > 3.3TeV). 
If a stop is significantly lighter than the other squarks the corresponding limits become 
rrif- > 950 GeV (m^ > 2 TeV). In both cases the gluino is expected to be either excluded 
or discovered up to a mass of nig ~ 1.4 TeV {rrig ~ 2.6 TeV). Intermediate cases may be 
estimated by interpolating the results according to the discussion in section 12.11 Already 
the data collected so far should allow one to place lower limits above a TeV on both 
rUq and rrig even in the most challenging scenario. The actual experimental searches 
might achieve even better sensitivities, since we chose quite conservative assumptions 
concerning background rejection and detector efficiencies. 

Staus stopped in the detectors could provide intriguing possibilities to test super- 
gravity or to gain insights into the SUSY breaking mechanism. We have computed the 
number of stopped stau events in the framework of the simplified models. Especially 
for very compressed spectra, the LHC provides a very good environment to measure the 
stau lifetime. 
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